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Abstract 

We calculated the decay widths of the l^i<2 and l^i<4 ss mesons and com- 
pared them to the measured properties of the .^(2220) (now known as the 
/4(2220)). Including previously neglected decay modes we found that the 
width of the state ss meson is much larger than previously believed mak- 
ing this explanation unlikely. On the other hand the predicted width of the 
^i*4 state, although broader than the observed width, is consistent within the 
uncertainties of the model. This interpretation predicts large partial widths 
to KK*{^^2) and iT* (892)ir* (892) final states which should be looked for. 
A second possibility that would account for the different properties of the 
^(2220) seen in different experiments is that two hadronic states exist at this 
mass. The first would be a broader ^F^ ss state which is seen in hadron beam 
experiments while the second would be a narrow state with high glue content 
seen in the gluon rich J radiative decay. Further experimental results are 
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needed to sort this out. 
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I. INTRODUCTION 



It is roughly a decade since the ^(2220), now known as the /4(2220), was discovered by 
the MARK III collaboration in J/ip radiative decays to K^K^ and KsKs final states [|l|. 
Its most interesting property, which attracted considerable attention, was its narrow width 
of roughly 30 MeV. Because the width was inconsistent with expectations for a conventional 
qq meson with such a large mass, the ^'s discovery led to speculation that it might be a Higgs 
boson 0, a bound state of coloured scalars 0], a four quark state 0,0, a AA bound state 
0, a hybrid or a glueball Despite the prevailing wisdom, the authors of Ref. |^,|^ 
argued that the properties of the ^(2220) could be consistent with those of a conventional 
meson: the L=3 ss meson with J^^ = 2++ or J^'-'' = 4++. 

In the original analysis of L=3 ss properties it was shown that of the qq states with the 
appropriate J^'-' quantum numbers only the and ss states of the first L=3 multiplet 
have masses consistent with the ^(2220) According to this analysis these two states 
were exceptional in that they have a limited number of available decay modes which are all 
relatively weak. However, the analysis was not exhaustive in that it did not calculate the 
decay widths to all possible final states. In particular it made the assumption, which we 
will see to be incorrect, that the decays to an L = 1 meson and a. K or rj were small on the 
basis of phase space arguments alone. 

To further complicate the discussion, more recent experiments have observed a hadronic 
state decaying to KK in different reactions and with different properties. The various 
experimental results relevant to the ^(2220) are summarized in Table |. The most recent 



measurement of the ^(2220) properties by the BES collaboration ||TT| indicates that its 
decays are approximately fiavour symmetric giving support to the glueball interpretation. 
At the same time, although the narrow ^ (2220) was not seen in J/ ip radiative decays by the 
DM2 experiment despite the fact that DM2 has slightly higher statistics, DM2 did observe 
a broader state decaying into KK [|10|. If all the experiments are taken at face value the 
overall picture is confused and contradictory. 
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In this paper we re-examine the nature of the ^(2220)//4(2220) meson and calculate the 
partial widths of the and ss states to all OZI-allowed 2-body final states allowed 
by phase space. To give a measure of the reliability of our analysis we calculate the widths 
using both the ^Pq decay model (often referred to as the quark-pair creation decay model) 
T^,|TB| and the flux-tube breaking decay model []T^. As an additional consistency check we 



calculated several partial widths using the pseudoscalar decay model [|T^]. Our goal is to 
shed some light on the nature of the ^ (2220) by comparing the quark model predictions for 
the hadronic widths to the various experimental results. 

The outline of the paper is as follows. In section |I| we briefly outline the models of 
hadron decays and the fitting of the parameters of the models. We relegate the details 



to the appendices. In section |T| we present the results of our calculations for the L = 
3 mesons and discuss our results. In the final section we attempt to make sense of the 
various contradictory experimental results and put forward our interpretation along with 
some suggested measurements which may clear up the situation. 

II. MODELS OF MESON PROPERTIES AND DECAYS 

The quark model has proven to be a useful tool to describe the properties of hadrons. 
The quark model has successfully described weak, electromagnetic, and strong couplings Q. 
In some cases we will use simplified meson wavefunctions which have been used elsewhere to 
describe hadronic decays Jl^ while in other cases we will use more complicated wavefunctions 
from a relativized quark model which includes one-gluon exchange and a linear confining 
potential [|18|. The strong decay analysis was performed using the QCD based flux-tube 



breaking model |T3- It has the attractive feature of describing decay rates to all possible 
final states in terms of just one fitted parameter. We also include results for the ^Pq model, 



often referred to as the quark-pair creation model pr5|Jl6|] , which is a limiting case of the flux- 



^See for example Ref. |18]. 
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tube breaking model and which greatly simplifies the calculations and gives similar results. 
As a final check we calculated some partial widths using the pseudoscalar emission model 
p!8| and confirmed that it also gave results similar to those of the fiux-tube breaking model. 



A. Decays by the ^Pq Model 



The ^Pq model |]T5|,|l6l is applicable to OZI-allowed strong decays of a meson into two 



other mesons, as well as the two-body strong decays of baryons and other hadrons. Meson 
decay occurs when a quark-antiquark pair is produced from the vacuum in a state suitable 
for quark rearrangement to occur, as in Fig. ||. The created pair will have the quantum 
numbers of the vacuum, ^Pq. There is one undetermined parameter 7 in the model - it 
represents the probability that a quark-antiquark pair will be created from the vacuum. 
The rest of the model is just the description of the overlap of the initial meson (A) and the 
created pair with the two final mesons (B,C), to calculate the probability that rearrangement 
(and hence decay) will occur. A brief description of the model is included in Appendix 
and the techniques by which the calculations were performed are discussed in Appendices 
g and g. 

B. Decays by the Flux- Tube Breaking Model 

In the fiux-tube picture a meson consists of a quark and antiquark connected by a tube 
of chromoelectric fiux, which is treated as a vibrating string. For mesons the string is in its 
vibrational ground state. Vibrational excitations of the string would correspond to a type 
of meson hybrid, particles whose existence have not yet been confirmed. 



The fiux-tube breaking decay model [O] is similar to the ^Pq model, but extends it 



by considering the actual dynamics of the fiux-tubes. This is done by including a factor 
representing the overlap of the fiux-tube of the initial meson with those of the two outgoing 
mesons. A brief review of the model is given in Appendix ^ and the techniques by which 
the calculations were performed are discussed in Appendices and 0. 



C. Fitting the Parameters of the Decay Models 



The point of these calculations is to obtain a reliable estimate of the and ss 
meson decay widths. To do so we considered several variations of the flux-tube breaking 
model. By seeing how much the results vary under the various assumptions we can estimate 
the reliability of the predictions. 

The first variation lies with the normalization of the mock meson wavefunctions and the 
phase space used to calculate the decay widths |T^. In the Appendices we have normahzed 
the mock meson wavefunctions relativistically to 2E and used relativistic phase space, which 
leads to a factor of EbEc/Ma in the final expression for the width in the centre of mass 
frame. We will refer to this as relativistic phase space/normalization (RPSN). However, 



there are arguments ||2^ that heavy quark effective theory fixes the assumptions in the 
mock meson prescription and suggests that the energy factor be replaced by MbMc/Ma, 
where the Mj are the calculated masses of the meson i in a spin-independent quark- ant iquark 
potential [|17|. (In other words Mi is given by the hyperfine averaged mass that is equal to 
the centre of gravity of the triplet and singlet masses of a multiplet of given L.) We will 
refer to this as the Kokoski-Isgur phase space/normalization (KIPSN). 

The second variation in our results is the choice of wavefunctions. We calculate decay 
widths for two cases. In the first we use simple harmonic oscillator (SHO) wavefunctions 
with a common oscillator parameter for all mesons. In the second case we use the wave- 



functions, calculated in a relativized quark model, of Ref. |I8| which we will label RQM. 
In all we looked at six cases: the ^Pq model using the SHO wavefunctions, the flux-tube 
breaking model again using the SHO wavefunctions, and the flux-tube breaking model us- 



ing the RQM wavefunctions of Ref. |T8[; in all three cases we used both choices of phase 
space/normalization. 

Some comments about the details of the calculations are in order. For the SHO wavefunc- 
tions, we took for the oscillator parameter (3 = 400 MeV which is the value used by Kokoski 



and Isgur ||T^. However, different quark models find different values of (3 so that there is 
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the question of the sensitivity of our results to (3. We will address this issue below. We used 
quark masses in the ratio m„ : rrid '■ mg = 3:3:5 — this differs from the calculations of 
Ref. IJT^, which ignored the strange-quark mass difference. In the RQM wavefunctions these 



parameters are already set — the values of (3 were found individually for each meson, and 
the quark masses were fitted: m„ = 220 MeV, = 220 MeV, and = 419 MeV. We have 
treated all mesons as narrow resonances, and have ignored mass differences between members 
of the same isospin multiplet 0. Masses were taken from the Review of Particle Properties 
1994 if the state was included in their Meson Summary Table 0. If it was not, then 



the masses predicted in Ref. [T^ were used. (This includes the masses of the 1^F2 and l^F^ 
ss mesons: 2240 MeV and 2200 MeV respectively.) Meson flavour wavefunctions were also 
taken from Ref. |]18[ - for the isoscalars we assumed ideal mixing (0nonstrangc = '^i.'^^ + 



Estrange = ss), cxccpt for the radial ground state pseudoscalars, where we assumed perfect 

mixing {jprj ■^^(0nonstrange Estrange) > 0r)' (0nonstrange ~l~ Estrange))- 

We fitted 7, the one undetermined parameter of the model, in a global least squares 
fit of 28 of the best known meson decays. (We minimized the quantity defined by = 
^^j^-pmodei _ -peip^2^^2^ whcrc (Jr, is the experimental error Q.) The experimental values for 
these decays and the fitted values for the six cases are listed in Table To give a more 
descriptive picture of the results we plotted in Fig. ^, on a logarithmic scale, the ratio of 
the fitted values to the experimental values. From Table |I| one can see that the results for 
the ^Pq and fiux-tube breaking models for the SHO wavefunctions are very similar^. We 



^The one exception was for the decay (/> K~^K where the charged and neutral kaon mass 
difference is significant to the phase space. 



^The one exception was the I^Pq ss state — see Table IV . 

^For the calculations in the flux-tube breaking model, a 1% error due to the numerical integration 
was added in quadrature with the experimental error. 
^The one exception to this is the S-wave decay -ftrQ(1430) Kn which seems particularly sensitive 
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therefore only plotted the ^Pq model results using the SHO wavefunctions and the flux-tube 
breaking model results for the RQM wavefunctions. A reference line is drawn in each case 
for p'^o'^ey-pexp _ ]^ guide the eye. Since all the partial widths are proportional to 7^, 
using a different fit strategy rescales 7. This is equivalent to simply shifting all points on the 
plot simultaneously making it easy to visualize any change in agreement for specific decays. 

The KIPSN gives a better overall fit to the data. Even so, certain decays, ^3(1780) — >■ 
Kp and /4(2050) — > uu for example, are fit much better using the RPSN. For both the 
RPSN and KIPSN one can see in Fig. I that a significant number of the decays differ from 
the experimental values by factors of two or more. Decays with two pseudoscalars in the final 
state tend to do better with the KIPSN but the KIPSN generally underestimates decays of 
high L mesons with vector mesons in the final states. On the other hand the RPSN tends 
to overestimate decays with two pseudoscalars in the final states. Similar observations can 
be made for the fiux-tube breaking model using the RQM wavefunctions. Having said all 
this we stress that these are only general observations and exceptions can be found to any 
of them in Table y. One must therefore be very careful not to take the predictions at face 
value but should try if possible to compare the predicted decay to a similar one that is 
experimentally well known. 

Finally, we consider the sensitivity of our results to p. In addition to the fits discussed 
above, we performed simultaneous fits of both 7 and (3 to the 28 decay widths for both the 
RPSN and the KIPSN. The resulting values of 7 and (3 are 13.4 and 481 MeV respectively 
for RPSN and 5.60 and 371 MeV respectively for KIPSN. In both cases the overall fits 
improved slightly, with some widths in better agreement and some in worse agreement with 
experiment when compared to the fits for f3 = 400 MeV. However, the fitted widths of the 
most relevant '^^4 decays improve slightly for RPSN but show mixed results for KIPSN. 
We also redid our fits of 7 to the decay widths for (3 = 350 GeV and (3 = 450 GeV. For 



to the model. 
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j3 = 350 MeV the overall fit improves slightly for KIPSN although the predicted /4(2050) 
decay widths are a little worse and the i^'4(2045) widths are a little better. For RPSN the 
overall fit is a little worse as are the decays. For (3 = 450 MeV the overall fit with KIPSN 
becomes a little worse as does the fitted ^F^ widths while for RPSN the overall fit and fitted 
^^4 widths become a little better. We conclude that while there is some sensitivity to /?, the 
results for modest changes in /3 (including the /? we obtain by fitting 7 and P simultaneously) 
are consistent with those for f3 = 400 MeV within the overall uncertainty we assign to our 
results. It should be stressed that it is not sufficient to simply change f3 but that a new 
value of 7 must be fitted to the experimental widths included in our fit. 

III. RESULTS FOR AND ^^4 ss MESON DECAYS 

Using the 7's obtained from our fit we calculated all kinematically allowed partial widths 
for the and ^^4 ss meson decays. The results are given in Tables ^ and [rV| . 
For the ^^4 state the main decay modes are: 

/4 ^ ir*(892)ir*(892), KK, KK*{892), (f)(j), KK;{im), KKi{im), r]r], r]r]' (1) 

For the KIPSN and the SHO wavefunctions the total width is 132 MeV with the ^Pq 
model. For this set of assumptions the KK, rji], and rjr]' modes are probably reasonably 
good estimates. However, the decay widths to KK* (892) and K* (892) K* (892) are likely to 
be larger than the predictions. On this basis it does not seem likely to us that the /4 width 
is less than the predicted total width by a factor of two or more, i.e. we do not expect it to 
be less than about 70 MeV. If anything, we would expect it to be larger than the predicted 
width, i.e. > 140 MeV. 

For the state we obtain results similar to the ^^4 state for the KK, KK* (892), 
and K* {892) K* {892) modes. However, the also has large partial widths to KKi{1270), 
/s:*(892)i^'i(1270), KK*{U30) and r//i(1510). In fact, KKi{1270) is the dominant decay 
mode. It is large in all variations of the calculation we give in Table IV. The most closely 
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related decay in our fit is the decay 7r2(1670) /2(1270)7r which is relatively large and is 
well reproduced by the KI normalization and SHO wavefunction case. The total width for 
this case is ~ 400 MeV Even if this width is overestimated by a factor of two, it would 
still be too large to identify with the ^(2220). 

Although this result appears surprising it has a straightforward explanation. Examining 
Table |V[ the lowest angular momentum final states in /2 decay are P-waves. All of these 
decays are relatively broad but the Ki{1270)K is the P-wave decay with the largest 

available phase space. In fact, one could almost order the P-wave decays using phase space 
alone. The analogous decay of the /4 is in an F-wave and therefore is subject to a larger 
angular momentum barrier. The lowest angular momentum partial wave for decays is a 
D-wave which although it has the largest partial width of all decays is still smaller than 
the P-wave decay. 

As another measure of the reliability of these predictions we calculated the widths of 
the ^^'4(2045) and /4(2050) mesons (the ^^4 i^'-like and non-strange isovector mesons, re- 
spectively). The results for all significant kinematically allowed final states are given for 
the ^Pq model using SHO wavefunctions in Tables |V| and |VI| respectively. The results are 
consistent with the general fit results given in Table ^ and Fig. |^. In general, the widths 
calculated using RPSN tend to be larger and those calculated using the KIPSN tend to 
be smaller. More specifically, decays to two pseudoscalar mesons using RPSN are gener- 
ally overestimated while the results calculated using KIPSN are in reasonable agreement 
with experiment. There is no pattern for the decays to two vector final states. The decay 
-^'4(2045) K*{892)p is greatly overestimated using RPSN but is in good agreement us- 
ing KIPSN. In contrast, the predicted decay /4(2050) ujuj agrees well using RPSN but 

^We note that the LASS collaboration has observed a i^'KlOSO) state with a large total width of 
373 lb 33 lb 60 MeV which could be associated with the strange meson partner of the ^^2(55) meson 

m 
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is greatly underestimated using KIPSN. The total widths tend to be overestimated using 
RPSN but are underestimated using KIPSN, both to varying degrees. The only conclu- 
sion we can draw from these results is that the total width probably lies between the two 
estimates but it is difficult to guess if it is closer to the lower or upper value. 

Finally, in Table VII we give the predicted total widths for the and ^F^^ ss states for 
the different values of (3 considered in the previous section. Although they vary considerably, 
by roughly a factor of 2 going from (3 = 350 MeV to /3 = 450 GeV (except for the r{^F2) with 
RPSN which varies by a factor of 3), these values are consistent within the large uncertainties 
we assign to our results. 

IV. DISCUSSION AND CONCLUSIONS 

The motivation for this paper was to re-examine the possibility that the ^(2220) is an 
L=3 ss meson. This question is especially timely given the recent BES measurements of 
a narrow resonance with a mass of 2.2 GeV seen in J/ip radiative decays. To do so we 
calculated all kinematically allowed hadronic decays of the and '^^4 ss states using 
several variations of the flux-tube breaking decay model. 

It appears very unlikely that the ^(2220) can be understood as the ^^2 ss state. All 
variations of our calculation indicate that the ss is rather broad, ^ 400 MeV. The 
dominant decay mode is the difficult to reconstruct KKi{1270) final state. Other final states 
with large branching ratios are /s:*(892)is:i(1270), KK*{892), K* (892) K* (892), KK*{U30), 
KK, and r//i(1510). 

It is more likely that the ^^4(55) state can be associated with the ^(2220). The calculated 
width is ~ 140 MeV but given the uncertainties of the models it is possible, although perhaps 
unlikely, that the width could be small enough to be compatible with the width reported 
in the Review of Particle Properties 1994 [^. In this scenario the largest decay modes 
are to K* (892) K* (892), KK, KK* (892), and (jxf). Since only the KK final state has been 
observed an important test of this interpretation would be the observation of some of these 
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other modes. 

There are, however, some problems with the ^^4(55) identification of the ^(2220). Fore- 
most is the flavour symmetric decay patterns recently measured by the BES collaboration 



Tl| . These results contradict the expectations for a conventional ss meson. Second is the 
wide range of measured widths for this state. Although the Review of Particle Properties 
1994 lists an average width of 38^ J| MeV the widths measured in hadron production exper- 
iments, LASS and E147, are larger while those measured in J/ip radiative decay tend to be 
narrow. The exception is the DM2 experiment which does not see, in J/ip radiative decay, 
a narrow state in J ftp ^ jKK but does observe a relatively broad state at this mass. 

To account for these contradictions we propose a second explanation of what is being 
observed in this mass region — that two different hadron states are observed, a narrow state 
produced in J/ ip radiative decay and a broader state produced in hadron beam experiments. 
The broader state would be identified with the ^F^^ss) state. The predicted width is consis- 
tent with the quark model predictions and the LASS collaboration shows evidence that its 
quantum numbers are J^'" = 4++. We would then identify the narrow hadron state observed 
in the gluon rich J/ip radiative decays as a glueball candidate predicted by lattice gauge 



theory results |]22|. Recent lattice results indicate that glueballs may be narrower than one 
might naively expect The scalar glueball width is expected to be less that 200 MeV 



and one might expect a higher angular momentum state to be even narrower. The narrow 
state is not seen in hadron beam production because it is narrow, is produced weakly in 
these experiments through intermediate gluons, and is hidden by the ss state. Conversely, 
the broader state is not seen in J/ ip radiative decays since this mode preferentially produces 
states with a high glue content. Crucial to this explanation is the experimental verification 
of the BES results on the flavour symmetric couplings of the state produced in J/ip radia- 
tive decay and the observation of other decay modes for the broader state in addition to 
the theoretical verification that the predicted tensor glueball is as narrow as the observed 
width. 

The ^(2220) has been a longstanding source of controversy. It is a dramatic reminder that 
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there still is much that we don't understand about hadron spectroscopy and demonstrates 
the need for further experimental results to better understand this subject and ultimately 
better understand non-Abelian gauge theories, of which QCD is but one example. 
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APPENDIX A: REVIEW OF THE ^Pq MODEL OF MESON DECAY 

We are looking at the meson decay A BC in the ^Pq model (Fig. [ip. Define the S 
matrix 

S = I -27^1 6{Ef - E,) T 

and then 

{f\T\i) = 5\Pf - Pi) M^^^A^^'B^Jc (Al) 
which gives, using relativistic phase space, the decay width in the centre of mass (CM) frame 

^ 2 P S 



TC 



\M^Ja^JbMJc\^. (A2) 



Here P is the magnitude of the momentum of either outgoing meson. Ma is the mass of meson 
A, \ Ja, Mj^) are the quantum numbers of the total angular momentum of A, s = 1/ {1 + 6bc) 
is a statistical factor which is needed if B and C are identical particles, and M'^^-'a^'^Jb^'^-'c 
is the decay amplitude. 

For the meson state we use a mock meson defined by 



\A{nA''-^'LA ){Pa)) = PEa E {LaMl^SaMsAJaMj^) 

X / dV^ ^uaLaMi^J^Va) xSms, 0a 



X 



\(1i{i^,Pa+Pa) q2{^PA - pa)) (A3) 
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The subscripts 1 and 2 refer to the the quark and antiquark of meson A, respectively; pi 
and nil are the momentum and mass of the quark. Note that the mock meson is normahzed 
relativistically to 2Ea S^{Pa — P'a)i but uses nonrelativistic spinors and CM coordinates 
{Pa = Pi + P2 is the momentum of the CM; pA = {m2Pi — 1TL1P2) / {mi + 1712) is the relative 
momentum), ua is the radial quantum number; \La,Ml^) and \Sa,Ms^) are the quantum 
numbers of the orbital angular momentum between the two quarks, and their total spin 
angular momentum, respectively; {LaMl^SaMs^\JaMj^) is a Clebsch-Gordan coefficient. 
Xs^Afg^, 0A aiid Ua are the appropriate factors for combining the quark spins, flavours 
and colours, respectively, and iPuaLaMlJ^Pa) is the relative wavefunction of the quarks in 
momentum space. 

For the transition operator we use 

T = -37 E {Iml-mm I dVa <^'p,5\n + Va) yT{^) xl\n € ^0' b\{p,) d\{p,) 

(A4) 

where 7 is the one undetermined parameter in the model ^ and yi^{p) = p''Yi"^{9p, (pp) is 
a solid harmonic that gives the momentum-space distribution of the created pair. Here the 
spins and relative orbital angular momentum of the created quark and antiquark (referred 
to by subscripts 3 and 4 respectively) are combined to give the pair the overall J^*" = 0+"'" 
quantum numbers (in the ^Pq state). 

Combining Eq. ^3| and gives for the amplitude in the CM frame (after doing the 
colour wavefunction overlap): 

M^j^Mj^Mj^^p^ = 7 ^SEaEbEc J2 (LaMlaSaMsJJaMj^) 

MlaMSaMLsMSb, 



^Ouv value of 7 is higher than that used by Kokoski and Isgur |[T^ by a factor of VSOvr due to 
different field theory conventions, constant factors in T, etc. The calculated values of the widths 
are, of course, unaffected. 
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X {LbMl^SbMs^ I JBMj^){LcMLaScMsc \ JcMj^){lm 1 -m|00) 



+ (-1) 



1+5^+5b+5c/j,32 ,14, ,12 ,34 



(A5) 



The two terms in the last factor correspond to the two possible diagrams in Fig. |I| - in the 
first diagram the quark in A ends up B; in the second it ends up in C. The momentum space 
integral I {P, mi, 1712,1713) is given by 

I{P, rm,m2, ms) = / d'prn.LsM^^/ + P) rncLcM^^/ + P) ^n^L^M, + P) y^iv) 

(A6) 



where we have taken P = Pb = —Pc- 



APPENDIX B: REVIEW OF THE FLUX-TUBE BREAKING MODEL OF 

MESON DECAY 

The flux-tube breaking model of meson decay extends the ^Pq model by considering the 
actual dynamics of the fiux-tubes. This is done by including a factor representing the overlap 
of the fiux-tube of the initial meson with those of the two outgoing mesons. Kokoski and 
Isgur have calculated this factor by treating the fiux-tubes as vibrating strings. They 
approximate the rather complicated result by replacing the undetermined parameter 7 in 
the ^Pq model with a function of the location of the created quark-antiquark pair, and a 
new undetermined parameter 70: 

Here b is the string tension (a value of 0.18 GeV^ is typically used) and Wmin is the shortest 
distance from the line segment connecting the original quark and antiquark to the location 
at which the new quark-antiquark pair is created from the vacuum (see Fig. |^): 

ifr>w|cos6'| 

w ■ = < 

mm I 

_|_ _ |cos6'| , if r < w |cos^| 
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To incorporate this into the ^Pq model, we first Fourier transform Eq. ^6| so that the 
integral is over position-space. We then pull the parameter 7 inside the integral, and replace 
it by the function of position 7(r, tv). The expression for the amplitude in the flux-tube model 
is then the same as that of Eq. ^ except that 7 is replaced by 70, and /(P, mi, 1712, m^) is 
replaced by 



mi, 1712, ms 



{2n 



dfdw ^PI^LbMlJ-^ - r) V 



tlcLcMl, 



{w — r) 



fj^=—2r 
™2 



g 2 mill 



where the ip^s are now the relative wavefunctions in position space. 



APPENDIX C: CONVERTING TO PARTIAL WAVE AMPLITUDES 

The decay amplitudes of the ^Pq and flux-tube breaking models derived in Appendices |^ 
and 0, M'^-^a^'^Jb^^Jc ^ are given for a particular basis of the final state: \9, (J), Mj^, Mj^) = 
\Q, Mjg, Mj^). Here 6 and are the spherical polar angles of the outgoing momentum of 
meson B in the CM frame. 

We would prefer to calculate amplitudes for particular partial waves, since they are what 
are measured experimentally: \J,M,S,L). Here \J,M) are the quantum numbers of the 
total angular momentum of the final state, \S,Ms) are the quantum numbers for the sum 
of the total angular momenta of B and C, and \L,Mi) are the quantum numbers for the 
orbital angular momentum between B and C. 

The formula for the decay width in terms of partial wave amplitudes is different from 
Eq. ^: 

S,L 

where 

-pSL ^ ^ I T\/rSL\2 

A MX ' ■ 
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is a partial wave amplitude, and V^^ is the partial width of that partial wave. 
We used two methods to convert our calculated amplitudes to the partial wave basis 
a recoupling calculation, and by use of the Jacob- Wick Formula. 



P5| 



1. Converting by a Recoupling Calculation 

The result of a recoupling calculation is 

(CI) 

Note that this can be done for any value of Mj^ ; alternatively, one could sum over Mj^ and 
divide by {2 J a + 1), on the right side. 



2. Converting with the Jacob- Wick Formula 

The Jacob- Wick formula relates the partial wave basis | J, M, S*, L) to the helicity basis 
|J, M, Ab,Ac), where and Ac are the hehcities of B and C, respectively. To use it we 
must first convert the basis that we calculate with to the helicity basis. This is done by first 
choosing P = Pb to lie along the positive z axis (in the CM frame still), so that A^ = Mj^ 
and Ac = —Mj^. Then one can use another expression that relates the helicity basis to the 
basis 1^2, A^, Ac)- 

The final result is 



,/47r(2L + 1) ^ 

M'\P)= E {L()S{Mj, + Mj^)\Jj,{Mj,+Mj,)){JbM.j,JcMj^\S{Mj, + Mj,)) 



X M'^^^^a =^Jb +Mjc Wjb ^^^Jc{Pz). 
Here Mj^ in the calculated amplitude is replaced by Mjg + Mj^ . 
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APPENDIX D: CALCULATIONAL TECHNIQUES 



The decay amplitudes in the ^Pq model were converted to partial wave amplitudes by 
means of a recoupling calculation. The whole expression for the amplitudes, including the 
integrals of Eqs. X6 and ^1|, was converted into a sum over angular momentum quantum 
numbers, using the techniques of Roberts and Silvestre-Brac (a result very similar 
to theirs was obtained). These techniques require that the radial portion of the meson 
wavefunctions be expressible in certain functional forms, which encompass simple harmonic 
oscillator wavefunctions. Our simple wavefunctions obviously meet these requirements, and 



since the detailed wavefunctions of Ref. |]18[ are expansions in terms of SHO wavefunctions, 
they do too. 

These expressions for the amplitudes were then computed symbolically using routines 
written for Mathematica [^. These routines are usable for any meson decay where the 
radial portion of the wavefunctions can be expanded in terms of SHO wavefunctions, and 
are limited only by the size of the symbolic problem that results, and the available computer 
resources. 

In the flux-tube breaking model there are two 3-dimensional integrations before convert- 
ing to partial wave amplitudes. The wish to be able to write general routines for any meson 
decay meant that only two of the six integrals could be done analytically; the remaining four 
must be done numerically. In order to minimize the numerical integration, the Jacob- Wick 
formula, rather than a recoupling calculation, was used to convert to partial wave amplitudes 
since no further integrals are involved. 

An integrand for each partial wave amplitude was prepared symbolically and converted 
to Fortran code using routines written for Mathematica, and then integrated numerically 
using either adaptive Monte Carlo (VEGAS or a combination of adaptive Gaussian 
quadrature routines. Again, these routines are usable for any meson decay where the radial 
portion of the wavefunctions can be expanded in terms of SHO wavefunctions, and are 
limited only by the size of the problem and available computer resources. 
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FIGURES 




FIG. 1. The two possible diagrams contributing to the meson decay A BC in the ^Pq model. 
In many cases only one of these diagrams will contribute. 
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FIG. 3. The position-space coordinates used in the flux-tube model. The cigar-shaped dashed 
line shows a possible surface of constant Wmin- 
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TABLES 



TABLE I. Summary of ^(2220) measurements. 



Experiment 


Mass 


Width 


Production 


Decays 




(MeV) 


(MeV) 






Mark III ^ 


2231 ± 8 


21 ±17 




BR{J/iP X BR{^ K+K-) 



= {4.2tl-l ± 0.8) X 10"^ 
BR{J/i^ ^ 70 X BR{C ^ KsKs) 

= {S.ltli ± 0.7) X 10-^ 
BR{J/ij 7^ X BR{(, vrvr) 

< 2 X 10-5 (90% C.L.) 
BRiJ/il^^jOxBRi^^pp) 

< 2 X 10-5 (90% C.L.) 

DM2 2230 26 J/ip^li BR{J/ip ^ 'j^) x BR{^ ^ K+K') 

< 2.3 X 10-5 (95% C.L.) 
BR{J/i; -> 7e) X BR{^ ^ i^^i^s) 

< 1.6 X 10-5 (95% C.L.) 

2197 ± 17 201 ± 51 J/ip jX BR{J/^P 7X) x BR{X KsKs) 

~ 1.5 X 10-^ 

BES<i 2233 ±5 19 ±11 J/^p^7$. BR{J/^p ^ j^) x BR{C ^ tt+tt-) 

= (5.61JI ± 1.4) X 10-5 

BR{J/ip 7^ X BR{£, pp) 

= (l.btoi ± 0-5) X 10-5 
BR{J/iP 7O X BR{^ K+K-) 

= (3.3l}| ± 1.1) X 10-5 
BR{J/i^ ^ 70 X BR{C ^ i^5i^5) 

= (2.7lJ:9 ± 1-0) X 10-5 

LASS^ 2209;j^± 10 60^^?^ R-p^ R-R+A 
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E147 f 2230 ±20 80 ± 30 vr^p KsKsn 

PS185 s 2231 30 pp KsKs BR{C pp) x BR{C KsKs) 

< 5.4 X 10"^ (3 S.D. J=4) 

^Ref. [|10| 

^Note that these values aren't measurements - they were assumed in order to set the BR Hmits. 
■^Ref. In I 



^Ref. PI 
fRef. 



gRef. Hi 
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TABLE II. Experimental and calculated widths (in MeV) of decays used in our global fit of 
the decay models' parameters. 



Decay r(experiment) ^Pq Flux- Tube Breaking Flux- Tube Breaking 

(SHO) (SHO) (RQM) 







RPSN KIPSN RPSN 


KIPSN 


RPSN 


KIPSN 


7 




9.73 


6.25 


16.0 


10.4 


20.5 


12.8 


P — ^ TTTT 


151.2 ± 1.2 


96 


148 


93 


148 


104 


152 


61 (1235) ^ COTT 


142 ±8 


176 


115 


155 


104 


306 


190 


02(1320) pn 


75.0 ±4.5 


65 


38 


67 


40 


84 


46 


02(1320) ^ KK 


5.2 ±0.9 


11 


8.0 


11 


8.5 


7.3 


5.0 


7r2(1670) ^ /2(1270)7r 


135 ± 11 


147 


116 


143 


117 


327 


246 


7r2(1670) pn 


74± 11 


232 


74 


226 


74 


323 


97 


7r2(1670) ^ K*{892)K + c.c. 


10.1 ± 3.4 


38 


17 


37 


17 


49 


21 


PsiimO) TTTT 


50.7 ±5.5 


116 


35 


122 


38 


68 


19 


P3(1690) OOTT 


34± 13 


36 


11 


39 


13 


45 


13 


/93(1690) ^ KK 


3.4 ±0.6 


9.2 


3.8 


9.7 


4.2 


4.2 


1.7 


/2(1270) ^ TTTT 


156.8 ± 3.2 


203 


109 


209 


116 


157 


80 


/2(1270) KK 


8.6 ±0.8 


7.2 


5.4 


7.4 


5.7 


5.0 


3.5 


(f) K+K- 


2.17 ±0.05 


2.37 


2.83 


2.28 


2.80 


2.30 


2.60 


/^(1525) ^ KK 


61 ±5 


117 


61 


118 


64 


98 


49 


K*{892) Ktt 


50.2 ±0.5 


36 


52 


34 


51 


38 


52 


K^{U3(}) Ktt 


267 ± 36 


163 


84 


117 


63 


875 


430 


1^1(1430) ^ Ktt 


48.9 ± 1.7 


108 


56 


112 


60 


88 


43 


i^|(1430) ^ i^*(892)7r 


24.8 ± 1.7 


27 


16 


27 


17 


31 


18 


ii:|(1430) ^ Kp 


8.7 ±0.8 


9.3 


4.9 


9.6 


5.2 


12 


5.8 


^^1(1430) ^ 


2.9 ±0.8 


2.6 


1.4 


2.6 


1.4 


3.2 


1.6 


i^|(1780) ^ i^p 


74± 10 


24 


7.7 


25 


8.4 


28 


8.7 
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Kg* (1780) ^ i^*(892)7r 


Ah ±7 


33 


11 


34 


12 


37 


12 


i^|(1780) ^ Kt: 


31.7 ±3.7 


87 


28 


92 


30 


54 


16 


ir|(2045) ^ Kit 


19.6 ±3.8 


55 


13 


59 


14 


28 


6.2 


K|(2045) ^ i^*(892)(/> 


2.8 ± 1.4 


3.2 


1.0 


3.3 


1.1 


4.7 


1.4 


/4(2050) ^ a;a; 


54± 13 


53 


11 


54 


11 


94 


18 


/4(2050) ^ TTTT 


35.4 ±3.8 


123 


25 


132 


28 


58 


11 


/4(2050) ^ KK 


1.4 ±0.7 


5.4 


1.6 


5.8 


1.7 


1.8 


0.5 
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TABLE III. Calculated partial decay widths (in MeV) for the ^F/^ ss state. We have calculated 
the widths of all kinematically allowed decays, but only show those partial widths that are > 1 MeV 
in at least one model. For this reason the total widths may not equal the sum of the partial widths 
shown. The subscripts on the decays refer to the S and L (see Appendix ^ of the given partial 
wave — the L is in spectroscopic notation (S, P, D, F, G, H). 



Decay 



3Pn 





(SHO) 



Flux- Tube Breaking Flux- Tube Breaking 



(SHO) 



(RQM) 





RPSN 


KIPSN 


RPSN 


KIPSN 


RPSN 


KIPSN 


fi ^ [KK]o,G 


118 


29 


125 


31 


62 


14 


fi ^ [KrK + c.c.]o,g'^ 


0.7 


0.4 


0.4 


0.2 


2.4 


1.2 


fi^[K*{892)K + c.c.]i,G 


107 


27 


115 


29 


112 


26 


fi^[K*{UW)K + c.c.]i,G'' 


1.7 


0.9 


0.8 


0.4 


5.0 


2.4 


fi ^ [Ki{mO)K + c.c.]iy 


6.4 


2.8 


7.0 


3.1 


10 


4.2 


fi ^ [Ki{mO)K + c.c^i^h" 


1.3 


0.6 


1.4 


0.6 


3.7 


1.5 


fi ^ [Ki{im)K + c.c.]iy 


14 


6.4 


15 


7.0 


29 


12 


fi^[K*^{U30)K + c.c.]2,F 


15 


7.0 


16 


7.7 


35 


15 


fi^[K*{892)K*{892)]o,G 


2.1 


0.5 


2.3 


0.6 


4.3 


1.0 


fi^[K*{892)K*{892)]2,D 


181 


44 


184 


46 


312 


72 


fi^[K*{892)K*{892)]2,G 


8.2 


2.0 


8.9 


2.2 


17 


3.9 


fi [w]o,G 


14 


3.5 


15 


3.9 


5.0 


1.2 


fi [v'v]o,G 


6.9 


1.7 


7.5 


1.9 


2.4 


0.6 


fi ^ [00]2,D 


20 


6.6 


21 


7.1 


31 


9.5 




498 


132 


522 


142 


633 


166 



"^Kr is our notation for the first radial excitation 

ii:*(i4io) 

K*(1680) 



'We used the following mixing p8| : < 



{2^ So) of the K. 

= 1.00(2^51) + 0.04(13Di) 

= -0.04(235i) + 1.00(13Di) 
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"We used the following mixing |21|: 



Ki(1270)+ = cos45°(liPi)+ + sin45°(l3Pi)+ 
Ki(1400)+ = -sin45°(liPi)+ + cos45°(l3Pi; 
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TABLE IV. Calculated partial decay widths (in MeV) for the ^i<2 ss state. We do not include 
a decay to /o(980)/o(980) because we question its assignment as a ^Pq QQ meson. At a more likely 
mass for the ^Pq ss meson, this decay is kinematically inaccessible. For other comments and notes, 

see Table g. 

Decay ^Po Flux- Tube Breaking Flux- Tube Breaking 

(SHO) (SHO) (RQM) 





RPSN 


KIPSN 


RPSN 


KIPSN 


RPSN 


KIPSN 


^ [KK]o,D 


51 


12 


47 


12 


101 


23 


^ [KrK + C.C.]o,D 


2.9 


1.5 


0.9 


0.5 


25 


12 




108 


26 


107 


26 


165 


38 


^ [K*{UW)K + c.c.]i,D 


2.6 


1.3 


0.6 


0.3 


4.0 


1.9 


^ [Ki{1270)K + c.c.]i,p 


445 


187 


449 


194 


1072 


426 


^ [Ki{mO)K + c.c.]i,F 


25 


11 


27 


12 


41 


16 


/^^ [Ki(1400)i? + c.c.]i,p 


14 


6.3 


15 


6.9 


29 


12 


f!, ^ [Ki{im)K + c.c.]i,F 


0.8 


0.4 


1.0 


0.4 


~ 


~ 


^ [K*2{U30)K + c.c.]2,P 


54 


24 


55 


25 


112 


47 


f^^[K*2{U30)K + c.c.]2,F 


9.6 


4.3 


10 


4.7 


22 


9.1 


f^^[K* (892) K*{892)]o,D 


24 


5.7 


24 


5.9 


39 


8.9 


f^^[K* (892) K*{892)]2,D 


14 


3.3 


14 


3.4 


23 


5.1 


/^^ [K*(892)K*(892)]2,G 


48 


12 


52 


13 


83 


19 


/^^ [Ki(1270)i?*(892)+c.c.]i,p 


99 


40 


102 


42 


209 


79 


/^^ [Ki(1270)i^*(892)+c.c.]i,p 


0.5 


0.2 


0.6 


0.2 


1.1 


0.4 


/^^ [Ki (1270)^* (892)+ c.c.]2,P 


33 


13 


34 


14 


70 


26 


^ [Ki(1270)i?*(892) +c.c.]2,F 


0.8 


0.3 


0.9 


0.4 


1.8 


0.7 


/2 ^ [w]o,D 


14 


3.3 


13 


3.2 


20 


4.4 




29 


7.0 


29 


7.2 


29 


6.6 


/^^ [/i(1510)7?]i,p 


45 


22 


46 


24 


92 


43 



30 



/2 ^ [/2(1525)77]2,P 14 6.9 14 7.3 29 14 

/2^['7V]o,D 6.6 1.6 6.7 1.7 4.9 1.1 

fi [#]o,D 3.9 1.2 3.9 1.3 5.5 1.6 

tl [(t>(t>h,D 2.2 0.7 2.3 0.7 3.1 0.9 

^ [#]2,G 1-0 0.3 1.0 0.3 1.1 0.3 

EiTj 1046 391 1058 406 2181 797 
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TABLE V. Calculated partial decay widths (in MeV) for the if|(2045) state. For comments 
and additional notes, see Table 



Decay Experiment ^Pq 

(SHO) 









RPSN 


KIPSN 


K|(2045) - 




19.6 ±3.8 


55 


13 


i^|(2045) - 


- [Kp]l,G 




19 


4.4 


K|(2045) - 


[K6i(1235)]i,F 




4.9 


2.2 


K|(2045) - 


[Kai(1260)]i,F 




1.3 


0.6 


K|(2045) - 


[Ka2(1320)]2,F 




2.2 


1.0 


K|(2045) - 


[K*(892)7r]i,G 




23 


5.5 


i^|(2045) - 


^ [K*{892)p]2,D 


> 18 lb 10 


76 


18 


Kl{2045) - 


^ [K*{892)p]2,G 




2.1 


0.5 


K|(2045) - 


[Ki(1270)7r]i,i. 




1.6 


0.7 


K|(2045) - 


[Ki(1400)7r]i,i. 




5.3 


2.6 


K|(2045) - 


[K2*(1430)^]2,F 




5.2 


2.6 


Kl{2045) - 






3.3 


0.9 


K|(2045) - 






6.0 


1.4 


i^|(2045) - 


[K4>]i,G 




1.1 


0.4 


K2{2045) - 


[K/ii(1170)]i,F 




2.9 


1.3 


K^(2045) - 


[K/2(1270)]2,F 




1.3 


0.6 


K|(2045) - 






4.9 


1.4 


ir|(2045) - 


[K*(892)w]2,D 




24 


5.7 


K|(2045) - 


[K*(892)(/.]2,D 


2.8 ± 1.4 


3.2 


1.0 






198 ± 30 


247 


65 



^This number is actually for the final state i^*(892)7r7r, and is the total for all partial waves. 
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TABLE VI. Calculated partial decay widths (in MeV) for the /4(2050) state. For comments 
and additional notes, see Table III, 



Decay 






Experiment 


RPSN 


(SHO) 

KIPSN 


/4(2050) - 






35.4 ±3.8 


123 


25 


/4(2050) - 


[7r7r(1300)]o,G 




3.9 


1.9 


/4(2050) - 


[7rai(1260)]i,F 




18 


7.5 


/4(2050) - 


[7ra2(1320)]2,F 




44 


19 


/4(2050) - 


[7r7r2(1670)]2,D 




2.1 


1.8 


/4(2050) - 


[pp]o,G 






1.9 


0.4 


/4(2050) - 


[pph,D 






159 


33 


/4(2050) - 


[pph,G 






7.3 


1.5 


/4(2050) - 


[vv]o,G 






3.2 


0.9 


/4(2050) - 


[vv']o,G 






1.0 


0.3 


/4(2050) - 


[77/2(1270)]2,F 




1.1 


0.5 


/4(2050) - 


■> [UJUJ]2,D 


> 


54 _|_ 1 Q a 


50 


10 


/4(2050) - 


■> [ww]2,G 






2.0 


0.4 


/4(2050) - 


[KK]o,G 




1.4 ±0.7 


5.4 


1.6 


/4(2050) - 


[KK*{892) + c.c.]i,G 




2.7 


0.8 


/4(2050) - 


[KKi{1270) +c.c.]i^F 




2.3 


1.2 


/4(2050) - 


[K*{892)K*{892)]2,D 




7.3 


2.1 








198 ± 30 


435 


109 



This number is the total for all partial waves. 
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TABLE VII. Calculated total decay widths (in MeV) for the and ss states for different 



values of /3. 



P (MeV) 


7 




r(3F2) 


r(3F4) 


RPSN 


350 


7.42 




590 


540 


400 


9.73 




1046 


498 


450 


12.0 




1549 


429 


48P 


13.4 




1841 


388 


KIPSN 


350 


5.16 




256 


170 


37P 


5.60 




309 


152 


400 


6.25 




391 


132 


450 


7.39 




534 


104 


^Prom the simultaneous fit of /? and 7. 
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